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ABSTRACT 

Context. X-ray data analysis have found that fairly complex structures at cluster centres are more common than expected. 
Many of these structures have similar morphologies, which exhibit spiral-like substructure. 

Aims. It is not yet well known how these structures formed or are maintained. Understanding the origin of these spiral-like 
features at the centre of some clusters is the major motivation behind this work. 

Methods. We analyse deep Chandra observations of 15 nearby galaxy clusters (0.01 < z < 0.06), and use X-ray temperature 
and substructure maps to detect small features at the cores of the clusters. 

Results. We detect spiral- like features at the centre of 7 clusters: A85, A426, A496, Hydra A cluster, Centaurus, Ophiuchus, 
and A4059. These patterns are similar to those found in numerical hydrodynamic simulations of cluster mergers with non-zero 
impact parameter. In some clusters of our sample, a strong radio source also occupies the inner region of the cluster, which 
indicates a possible connection between the two. Our investigation implies that these spiral-like structures may be caused by 
off-axis minor mergers. Since these features occur in regions of high density, they may confine radio emission from the central 
galaxy producing, in some cases, unusual radio morphology. 
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1. Introduction 

Previous studies have shown that many clusters con- 
tinue to be in the process of forming, groups or individ- 
ual galaxies being accr eated from the o uter large scale 
structure filaments (e.g. JPurret et al.ll200 3). Temperature 
maps of the X-ray emitting gas have shown that even 
clusters with apparently relaxed structures in X-ray can 
have perturbed temperature ma ps, which are indicative 
of recent ongoing mergers (e.g., Finoguenov et al.l l2005l : 
Durret fc Lima Netol 20081 ). The detection of these small 
structures is one of the main achievements of the present 
era of X-ray telescopes such as Chandra and XMM- 
Newton. However, the presence of cooling cores suggest 
that recent merger (s) have not had enough time to de- 
stroy the cooling core. 

More than 70% of cooling-flow clusters c ontain a cen- 
tral cD ga laxy that is a radio source fe.g.. lBurnslll990l : 
Eilekl [200! ) . Chandra observations of these systems found 
that the core of many clusters exhibit morphological com- 
plexities that are probably produced by the interaction 
between the intracluster medium (ICM) and the central 
radio galaxy. 

The most important connection between the X-ray 
ICM and the non-thermal relativistic electrons that pro- 



duce radio emission is dynamical. In the case of very en- 
ergetic jets, the ICM is heated and compressed because 



of su personic shocks (iHeinz et al ] Il998l : [Reynolds et al 



20011). However, there are weaker jets that do not pro- 



vide an efficient shock heating, but instead a bubble of 
low-density gas, which can be identified as a depression in 
the X-ray surface brightness, that rises because of buoy- 



ancy (e.g. JChurazov et al. ] l200ll : iBriiggen fc Kaiserll2002l : 



McNamara et al.ll2005] ). 



On the other hand, the ICM can also affect the ra- 
dio emission from radio galaxies by placing upper limits 
on the radio emission. The unusual radio morphology of 
central cluster galaxies that fill X-ray caviti es of low den- 
sity has been reported previously (Boehrin ger et al .11 19931: 



Tavlor et al][l994l : iFabian et al.ll200d I2OO2I : iTavlor et al 



20021 ). The most extreme case of radio distortion is the 



amorphous radio g alaxy P KS 07 45-1 91 at the core of the 
Centaurus cluster (| Tavlor et al.lll994h . 

Another interesting fea ture is the spiral- l ike feature at 
the core of Perseus (Churazov et al. 2003 ; Fabian et al.l 
2000l : lTavlor et al.ll2002l) A22Q4 (jSanders fc Fab ian 2008), 



and A2029 (jClarke et al.ll2004l ). This spiral pattern was 
also found in hydrodynamical simulati ons of merging clus- 
ters with non-zero impact parameters ([Gomez et al.ll2002 
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Ascasibar fc MarkevitchI 2QQ6I ). In these simulations, the 
infalling sub-cluster causes a disturbance in the mass peak 
as it passes through the cluster centre. Thus, the central 
cool gas of the main cluster acquires angular momentum, 
producing its spiral pattern. 

To understand more clearly the ICM-radio interaction 
and to characterize the spiral feature, we present a sample 
of 15 nearby clusters, among which 7 exhibit this feature in 
temperature and substructure maps, and a large cooling- 
flow and an embedded radio galaxy at the centre. The 
paper is organized as follows: in Sect. [21 we describe the 
data sample and the Chandra data analysis; in Sect.[3l we 
explain the image procedures used to detect the spiral-like 
pattern; in Sect. [H we present our results and notes for 
individual clusters; in Sect. [3 we discuss a possible sce- 
nario to explain the spiral-like feature formation and the 
ICM-radio connection; in Sect. [6l we present a synthetic 
image produced to reproduce the surface brightness of the 
Perseus cluster, and in Sect.[7J we summarize our findings. 

All distance-dependent quantities are derived by as- 
suming the Hubble constant Hq = 70 km s~^ Mpc~^, 
= 0.3 and, 1^a= 0.7. 



2. Observation and data reduction 

2.1. Cluster selection 

The objects in our sample were selected from a set of 
clusters within the redshift range 0.01 < z < 0.06 with 
Chandra public data. We imposed that the data for each 
clusters had a long exposure time, that is, of at least 35 ks. 
Our final sample contained 15 clusters, whose properties 
are presented in Tab. [H 



2.2. Chandra data reduction 

To obtain calibrated images without artifacts, that are ad- 
equate for achieving our goals, it is necessary to follow a 
series of procedures, otherwise we would have contamina- 
tion, which can be detected as substructures. 

We used the package CIAO 3.4. A level 2 event file 
was initially generated from a level 1 event file, using the 
standard procedur^il and the latest calibration of CALDB 
3.3.0. Periods with flares, which are spurious high counts 
caused by protons accelerated by the Sun, were excluded 
using the Ic-clean script. At this point, a rebinned image 
with pixels corresponding to 16 raw physical pixels (4x4) 
was created from the new level 2 event file, in the energy 
band from 0.3 to 7.0 keV. We then produced exposure 
maps and used them to obtain flat images from which 
the source points were removed by filling circles around 
each source by randomly sampling the same distribution 
as found in a circular region close to the source. Finally, 
we fitted a 2D analytical surface brightness model. 



http : / /cxc .harvard .edu /ciao3 .4/threads /createL2 / 



3. Imaging 

Since the bremsstrahlung emissivity depends on both the 
temperature and surface brightness, irregularities can be 
detected either in terms of temperature or density discon- 
tinuities. Since e ex T^/^ n^, these substructures are more 
significant in the surface brightness maps. Thus, to resolve 
the degeneracy, we used substructure (which depends on 
the surface brightness distribution, ?) and temperature 
maps to analyse in detail the inner parts of the clusters in 
our sample, of significant structure in the inner region of 
several clusters. To investigate this further, we compared 
these images with the temperature and the substructure 
maps. 

3.1. X-ray substructures 

We developed a new method to quantify substructures in 
clusters of galaxies using 2D surface brightness fits. This 
method was tested for 47 galaxy clusters observed with 
Chandra (?). 

The method is based on analysing the number and in- 
tensity of substructures detected by a friends- of -friends 
(FOF) algorithm. This analysis is performed on a resid- 
ual image, which is produced by fitting a bidimensional 
analytical model (/3-model or Sersic) with elliptical sym- 
metry (using Sherpa/CIAO, which fits a model by min- 
imizing the x^) to the cluster X-ray image. All surface 
brightness model parameters are free to vary, including 
the elipticity, position angle, and centre. In particular, to 
identify substructures, we used a threshold in the residual 
image to separate the pixels that had counts statistically 
significantly above or below the 2D fitted model surface 
brightness at the corresponding position. The threshold 
was defined to be 3 times the local X-ray background vari- 
ance. Once the pixels were separated, we identified these 
regions as substructures, using the FOF algorithm to link 
together the pixels above the threshold. This is a robust 
method if applied to high signal to noise ratio observa- 
tions, which is the case, since the clusters in our sample 
are rich objects at very low redshifts {z < 0.06) and were 
observed for more than 35ks. 



3.2. Temperature maps 

Ideally, one would prefer to compute a temperature 
map using: spatially resolved spectroscopy as we did in 
iLagana et al.l ( 2008 ). This, however, is very timing con- 
suming. We opt for another approach, namely to compute 
the hardness ratio maps and then converte them to tem- 
perature map. In this way, we are able to detect signifi- 
cant small-scale features without resorting to a computer 
intensive method. 

The hardness ratio (HR) maps, or color maps, are 
defined as the ratio of the fluxes in two different bands 
and can be interpreted as a projected temperature map. 
They can be defined as [smoothed hard-energy band im- 
age - smoothed soft-energy band image] / [smoothed hard- 
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Table 1: Observational information 



Cluster 


RA 
(J2000) 


Dec 
(J2000) 


z 


Detector 


Exposure time 
(ks) 


Soft band* 
(keV) 


Hard band* 
(keV) 


-j- 

njj 

(10^° cm-2) 


Abell 85 


00:41:37.80 


-09:20:33.0 


0.055 


ACIS-I 


38.91 


0.3-1.4 


1.4-6.0 


2.78 


Abell 426 


03:19:48.20 


+41:30:42.2 


0.018 


ACIS-S 


98.20 


0.3-1.4 


1.4-6.0 


13.6 


Abell 496 


04:33:37.00 


-13:14:17.0 


0.033 


ACIS-S 


76.08 


0.3-1.3 


1.3-6.0 


3.78 


Abell 3376 


06:02:10.00 


-39:57:21.0 


0.046 


ACIS-I 


44.85 


0.3-1.5 


1.5-6.0 


4.58 


Abell 754 


09:09:09.00 


-09:39:39.0 


0.054 


ACIS-I 


44.77 


0.3-1.5 


1.5-6.0 


4.82 


Hydra A 


09:18:05.70 


-12:05:42.5 


0.013 


ACIS-S 


100.13 


0.3-1.2 


1.2-6.0 


4.60 


Centaurus 


12:48:48.90 


-41:18:44.4 


0.011 


ACIS-S 


90.19 


0.3-1.2 


1.2-6.0 


8.56 


Abell 1644 


12:57:33.00 


-17:20:28.0 


0.047 


ACIS-I 


52.17 


0.3-1.6 


1.6-6.0 


4.11 


Abell 1991 


14:54:31.50 


+ 18:38:32.0 


0.059 


ACIS-S 


38.81 


0.3-1.1 


1.1-6.0 


2.46 


Abell 2052 


15:16:44.50 


+07:01:16.6 


0.035 


ACIS-S 


128.63 


0.3-1.2 


1.2 6.0 


2.71 


Abell 2107 


15:39:39.00 


+21:46:58.0 


0.041 


ACIS-I 


36.04 


0.3-1.5 


1.5-6.0 


4.45 


Ophiuchus 


17:12:27.80 


-23:22:11.5 


0.028 


ACIS-S 


51.18 


1.3-4.0 


4.0-7.0 


40.0 


Abell 3667 


20:13:07.25 


-56:53:24.0 


0.056 


ACIS-I 


105.01 


0.3-1.6 


1.6-6.0 


4.44 


Abell 4059 


23:57:00.70 


-34:45:33.0 


0.047 


ACIS-S 


93.34 


0.3-1.2 


1.2-6.0 


1.21 


Abell 2589 


23:23:57.40 


+ 16:46:39.0 


0.041 


ACIS-S 


54.13 


0.3-1.3 


1.3-6.0 


3.15 



* Soft and hard bands used to construct the hardness r atio. 

J Hydrogen column densities adopted to construct HR ( Dickey Lockmanlll990h . 



energy band image + smoothed soft-energy band image]. 
These bands are specified for each cluster in Table [1] and 
were chosen such that the two bands had almost identical 
net counts. 

We used the csmooth tool in Chandra Interactive 
Analysis of Observations (CIAO), an adaptive smoothing 
algorithm, to bin the data using bins of large angular size 
in order to examine possible structures in the image. In 
adaptive binning, the size of the smoothing kernel changes 
over the original image to create a constant S/N per pixel 
in the output image. The hard and soft energy band im- 
ages were smoothed identically, with the same adaptive 
kernel, avoiding spurious artifacts and yielding meaning- 
ful hardness ratios. 

To convert the HR maps into temperature maps, we 
required a simple algorithm to compare the relative fluxes 
in different X-ray bands with a theoretical plasma emis- 
sion model. Thus, based on the assumption of a single 
temperature, the HR-kT relation was estimated for each 
cl uster, and an examp le is shown in Fig. [T] for A496. 



Watanabe et al.l (|200ll ) also used this procedure to con- 
struct the temperature map of Ophiuchus using ASCA 
data. To do so, we took into account the weighted re- 
ponse matrices (RMFs), affective area files (ARFs) gener- 
ated with CIAO 3.4., exposure maps of each observation, 
the hydrogen column densities of each cluster (as specified 
in Table [2]), and adopted a fixed metallicity value of 0.35 
for all clusters. In Fig.[Tl we show that the HR conversion 
into plasma temperature exhibits little variation assuming 
different metallicity values. 




-0.6 -0.4 

Hardness ratio 

Fig. 1: Hardness ratio conversion into plasma temperature, 
using Abell 496 as an example. Top: HR-to-temperature 
conversion assuming different metallicity values (in solar 
units) when all other parameters remained fixed. Bottom: 
HR-to-temperature conversion assuming different hydro- 
gen column densities (in 10^^ cm~^). While a metallicity 
error has an insignificant effect, nn is important for an 
accurate determination of kT (in particular at high tem- 
peratures). 



4. Results 

In FiglSl we present the X-ray images, temperature and 
substructure maps of the 15 clusters in our sample. To in- 
vestigate the properties of these clusters, we also present in 



Tab |2] the mean X-ray temperature, the cooling- flow time, 
and the X-ray luminosity obtained from the literature. 

Although not clearly evident in all temperature maps, 
a spiral-like structure in the inner part of 7 clusters (A85, 
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Spiral-like structure at 



the centre of nearby clusters of galaxies 





Fig. 2: Cont. Left panels: X-ray image; middle panel: temperature maps; and right panels: substructure maps for A754, 
Hydra A cluster, Centaurus cluster, and A1644.The color-bar indicates the temperature in keV. 





Fig. 2: Cont. Left panels: X-ray image; middle panel: temperature maps; and right panels: substructure maps for 
A1991, A2052, A2107, and Ophiuchus cluster. The color-bar indicates the temperature in keV. 
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A426, A496, Hydra A cluster, Centaurus, Ophiuchus, and 
A4059) is detected in the substrucuture maps. These clus- 
ters are those of the highest cooling- flow rate and most are 
known to have a strong radio galaxy at the centre. 

Since the variations in temperature and density in- 
side shocks are positive correlated, that these patterns 
are also regions of low temperature and high density ar- 
gues against strong shocks. These substructures instead, 
resem ble the edges of cold fronts found in several c lusters 
(e.g., iMarkevitch fc Vikhlininl l2QQ7l : IPupke et all [2007). 
Cold fronts can be directly related to cluster merging 
activities, a nd by comparing our results with numerical 
simulations (|G6mez et al.l I2OO2I : lAscasibar fc MarkevitchI 
2006h one can see that the spiral-like features found here 



are morphologically similar to those found in their simula- 
tions. This indicates that their nature is probably dynam- 
ical and a possible scenario for its formation is discussed 
in Sect. \5\ 

Here, we provide additional information about each 
cluster in the sample, by highlighting key aspects of pre- 
vious studies. 

4.L Abell 85 

This i s a well-known sub clu mp cluster that has a cold- 
front ( Kempner et al.l [2002h and is a merger candidate. 
The optical properties of this cluster are indicative of 
evidence for galaxy star formation in the filament be- 
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Table 2: Cluster properties 



K. Arnaud^ 



Chen et al. (2007) f 



Ikebe et al. (2002)"f 



Edge et al. (1992) 



Cluster 



kT 
(keV) 



M 

(Mq yy- 



kT M -£>x[0-l - 2.4keV] 

(keV) (Mq yr-1) {Yv'^ lO^^erg/s) 



6-10-0.20 

3.80+°;- 



3 03 + 0-°^ 

7+0.60 
-0.60 
-, + 0.30 
-0.30 
3 70 + 2-20 



kT 
(keV) 



:[0.1 - 2.4keV] 
-2 1044erg/s) 



^ [2 - lOkeV] 
-2 1044erg/s) 



M 

(Mq yr-1) 



Abell 85 
Perseus (A426) 

Abell 496 
Abell 3376 

Abell 754 

Hydra A 
Centaurus 
Abell 1644 
Abell 1991 
Abell 2052 
Abell 2107 
Ophiuchus 
Abell 3667 
Abell 4059 
Abell 2589 



6.4 
5.7 



2.3 
5.1 
4.6 
3.7 
6.5 



128.5 
94.9 
100.7 



31.3 
19.3 
71.0 
30.8 
0.0 



7! 
4.10^ 



13.54 ± 0.15 

5.28 ± 0.07 
3.02 ± 0.07 
5.56 ± 0.15 
8.18 ± 0.06 
1.67 ± 0.06 
5.49 ± 0.48 

3.32 ± 0.06 



13.27 ± 0.15 
3.92 ± 0.08 
2.62 ± 0.06 



-0.23 
2+0.06 
-0.06 
+0.10 
-0.10 
.^ + 0.39 
-0.38 
9 00 + 0-"^^ 



4.59 
4.43^ 



3.12 



+0.10 
0.09 



10.25 
6.28 
3.94 
3.38 



+0.30 
0.36 
+ 0.27 
0.26 
+ 0.15 
0.15 
+ 0.13 
0.13 



1.81 ± 0.02 
2.93 ± 0.04 
0.69 ± 0.01 
0.74 ± 0.02 
4.0 ± 0.10 



0.73 ± 0.06 

0.43 ± 0.01 

2.2 ± 0.07 
1.78 ± 0.02 
0.52 ± 0.01 
3.44 ± 0.08 



10.56 
15.40 
3.56 

14.70 
3.89 
0.81 
3.33 

I. 97 

21.28 

II. 59 
2.65 



236 
183 
112 

24 

315 
18 
19 

90 

75 


124 



t All 
Ho = ' 



alues were converted to a common 
3 km s^-'^Mpc"-'- . 



cosmology with a Hubble constant given as 



tween the main cluster and the subclump (|BQue et al 
2QQ8I ), which is consistent with the previous hypothesis 
that the X-ray filament in Abell 85 is a gravitationally 
bound s tructure consisting of groups falling into th e main 
cluster (jLima Neto et alJl2QQlHDurret et alJl2QQ3h . 

We clearly see a spiral-like structure in the substruc- 
ture map. However, no corresponding structure appears 
in the temperature map. 

4.2. Abell 426 - Perseus 



Churazov et al.l (|2QQ3l ) presented XMM-Newton data in 
which a spiral-like structure was detected in temperature 
and surface brightness maps. These authors concluded 
that this feature is possibly a contact discontinuity sep- 
arating the main cluster gas from the gas of the infalling 
subcluster. A chain of galaxies is also associated with 
this region, probably tracing the filament along which the 
merger started. 

Since it has a strong central radio galaxy, NGC 1275, 
Boehringer et al. I (|i993) argued that the thermal plasma 
was displaced by the inner parts of the radio lobe. 
Fabian et al. ( 2QQ2I ) examined in detail the interaction be- 
tween the radio source 3C84 and the surrounding medium, 
arguing that the inner lobes are currently expanding sub- 
sonically. 

From our results, one can see that this cluster is one 
with the most prominent spiral-arm patterns in both tem- 
perature and substructure maps. The su bstructure found 
here is almost identical to that detected bv lChurazov et af 
( 20031 ) (see their Fig. 7). 



4.3. Abell 496 



This cluster hosts a prominent cold- front ( Dupke fc White! 
2003h . which is a sign of a recent minor merger, and 
although its resport has not received much attention. 



a very proni i nent spiral-like pattern was detected by 
iLagana et all (|2008h in Xym-Newton data. 

In Fig. O a spiral-like feature is evident in both the 
substructure and temperature map. The pattern that ap- 
pears in the temperature map is very similar in form to 
that of A85. 

4.4. Abell 3376 

This is a nearby ongoing merger clust er that was bserved 
by both Suzaku (iKawano et al. I l2008h and XMM- Newton 
(Bagchi e t al.ll2006[ ). A non-thermal radio-emitting struc- 
ture at the outskirts of this cluster was identified by the 
latter authors. They suggested that this structure prob- 
ably traces the elusive shocks of cosmological large-scale 
matter flows. 

4.5. Abell 754 

A754 was one of the first clusters for which an X-ray tem- 
perature map was derived from ROSAT (|Henrv fc Briel 
1995) and ASCA (jHenriksen fc MarkevitchI Il996h d ata. 



19961 ) 



in which the gas velocitv map ([Evrard et al 
dicated a plume-like pattern. A Chandra data analy- 
sis inferr a more complex merger g eometry, possibly 



involving more than two subcl usters (|Markevitch et al 



2003h. From X MM-ATe^^on data. lKassim et al.1 (|200l ') and 
Bacchi et al. confirmed that this cluster is a major- 

merger cluster and also detected a plume-like feature at 
its centre. There is a diffuse radio source associated with 



this cluster (jHenrv et al.ll2004[ ). 



In this work, a large substructure is detected but does 
not exhibit any particular shape. 

4.6. A780 - Hydra A cluster 

iFitchett fc MerrittI ( 19881 ) examined the dynamics of the 
Hydra cluster find that the velocity distribution of galaxies 
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close the cluster centre is very flat. Although they did not 
propose an infalling group scenario, they did consider the 
importance of substructures in this cluster, 

Chandra observations of the Hydra cluster detect a 
feature in the X-ray surface brightness surrounding radio 
lobes o f the act ive galactic nucleus (AGN) at the cluster 
centre (jNulsen et al. 2005) . These authors proposed that 
the shock front is driven by the expanding radio lobes in 
addition to the "new struct ure on small e r sca les" is asso- 
ciated with the radio source. IWise et aJ] (|2007 ) confirmed 
that the complex system found in this cluster was created 
by a continuous outflow or a series of bursts in the nucleus 
of the central galaxy. In this particular case, a merger has 
not been mentioned in the literature as a possible mean 
of forming this substructure. 

In this work, we were able to detect a truncated spiral- 
like feature in the substructure map that has a correspond- 
ing structure of relatively low temperature. 

4.7. Centaur us 

To explain the parti cular morphology of this cluster, 
Churazov et al.l ( 1999h analysed ASCA imaging and spec- 
tral results, and proposed that there is a subcluster cen- 
tred on NGC 4709 that is merging with the main cluster 
centred on NGC 4696. 

Sanders fc Fabian I |20Q6l) performed a detailed analy- 



sis of the core of the Centaurus cluster using a Chandra 
deep X-ray observation and XMM-TVewton data. To ex- 
plain the high metallicity present in the core of this clus- 
ter, the authors proposed a model in which the inner core 
of the Centaurus cluster has not experienced a major dis- 
ruption within the past 8 Gyr, or even longer. In their 
temperature map, the plume-like feature is evident, this 
finding then attracted little attention. They also reported 
the discovery of ripple-like X-ray surface brightness oscil- 
lations in the core of this cluster of galaxies, which they 
asserted were indicative of sound waves generated by the 
repeated inflation of central radio bubbles. 

From Fig. [21 one can see that this cluster has one of 
the most prominent spiral-like features detected in a sub- 
structure map. The very inner region of this enormous 
structure is also associated with a region of relatively low 
temperature in the temperature map. 

4.8. Abell 1644 

This is a complex merging system that was observed with 
XMM- Newto n, and found to cons ist of a main cluster and 
a subcluster ( Reiprich et al. 20041 ). These authors also de- 
tected a trail of cool, metal-rich gas closer to the subclus- 
ter. The combination of X-ray, optical, and radio data re- 
sults imply that the subcluster has passed by the main 
cluster off- axis. Although temperature maps were pre- 
sented in their work, the plume-like feature was not re- 
solved by XMM-Newton. 



In our analysis, a discontinuous spiral-like feature is 
present in the substructure map. This pattern is not as 
clearly evident as in Centaurus, Perseus, and A85 but it 
does resemble these structures. Unfortunately, only the 
inner part of the substructure has a corresponding feature 
in the temperature map. 

4.9. Abell 1991 



Sharma et~all (|2004l ) analysed XMM-TVewton data and 
detected an asymmetric surface brightness distribution 
with respect to the central galaxy. These authors de- 
tected bright knots of soft X-ray emission embedded in 
a cometary structure located north of the optical centre 
of the cD galaxy. According to these authors, the knots 
have no obvious association with the radio source. 

4.10. Abell 2052 



Zhao et al.l (|l993l ) reported a VLA detection of a radio 



Blanton et al.1 i2l 



sourc e at the core of this cluster. ^ 

carried out a Chandra X-ray data analysis of the 
large-scale properties of the cluster as well as the central 
region (which contains a radio source) presenting temper- 
ature and abundance profiles. 

Two clear substructures are present in Figl2]but they 
do not resemble to the other spiral-like feature presented 
in this work. 

4.11. Abell 2107 

iFujita et al.l (j2006 l) presented an analysis of Chandra ob- 
servation of this cluster. Since the cD galaxy has a large 
peculiar velocity and the ICM in the central region has 
an irregular structure, these authors concluded that this 
cluster is undergoing a merger. 

From our results, A2107 contains the smallest sub- 
structure and has almost symmetric X-ray contours pro- 
viding no evidence of a recent merger. 

4.12. Ophiuchus 

Ophiuchus is one of the hottest X-ra y cluster and is lo - 



cated 12 deg from the Galactic centr e ()Fu] it a et al.ll2008r) . 
Using ASCA data. iMatsuzawa et al.l ( 19961 ) acquir images 



in the energy range of 0.7-10 keV that indicated that the 
peak of the X -ray surface brightness is coincident with 
the cD galaxy. IWatanabe et al. (|200lh presented temper- 
ature and metallicity profiles obtained with ASCA data, 
which contained no evidence of a plume-like feature. In 
their analysis, the ROSAT archival data inferred an on- 
going merger, but a Chandra image indicates that this 
cluster has a cool-core that remains undisturbed. Thus, 
if the re is an ongoing; merg: er, it must be a minor merger 
fe.g.. lWatanabe et alJliool ). 

From the results presented in FiglSl Ophiuchus is the 
cluster that has the most prominent spiral-like pattern in 
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its substructure niap. However, no significant feature was 
detected in its temperature map. 

4.13. Abell 3667 

This cluster has a double radio halo outside its central re- 



gion (iRottgering et al. 1997; Ma rkevitch et al 1999). The 
galaxy distribution is bimodal (jProust et al.l Il988l ). the 
main component being located around the cD galaxy. 
A multiwavelength study confirmed that it is probably 
a merger cluster b ased on Chandra and optical data 
(jOwers et al.ll200^al ). Its X-ray image shows a clear cold 
front. 

4.14. Abell 4059 



Heinz et al.l (|2002l ) affirmed that there are clear signs of in- 
teraction between the radio galaxy and the ICM. However, 
these authors did not find any one-to-one spatial corre- 
spondence between the radio lobes and the X -ray cavities 
( Heinz et alll2QQ2[ ). iHuang SarazinI (|l998h argued that 
it is likely that the radio plasma from the cD galaxy may 
have displaced the X-ray gas and created the X-ray cavi- 
ties. 

On the other hand. lchoi et al.l ( 20041 ) analysed various 
scenarios concluding that the most probable is that of a 
compressio n front associated with ICM bulk motion. 

Reynolds_eL^lJ (2008) analysed a Chandra deep ob- 
servation, detect a subtle discontinuity in the gradient 
of the surface brightness that is approximately semi- 
circular in form and centred on the radio galaxy. T he 
Tevious C handra observatio ns reported by iHeinz et al.l 



and I Choi et al 



(|2004l ) were of insufficient signal- 
to-noise ratio to detect t his feature. Al though it can- 
not be rigorously proven, iRevnolds et al.l (|2008l ) inter- 
preted this feature as a weak shock caused by the radio 
galaxy activity. Similar shocks are seen in all hydro dy- 



namic simulations of jet /ICM i nteractions (|Heinz et al 
I2OO2I : IVernaleo R.evnQidsll2006l . e.g.). 



In this work, we have shown that this cluster also has 
a spiral-like substructure. The inner part of this pattern 
is related to a region of lower temperature, as it can be 
seen in the temperature map of FiglSl 

4.15. Abell 2589 

The hot gas in the core region of A2589 is undisturbed by 
interactions with a central radio source ( Buote fc LewisI 
2004h . Except for a 16 kpc shift in its X-ray cen- 
tre, A2589 has a rema rkably symmetrical X-ray image 
( Zappacosta et al]l2006l ). 



5. Discussion 

The literature contains detailed studies of large numbers 
of clusters with high cooling rates and a powerful radio 
galax y, e.g.. Hy dra A (McNamara et aDbOOd: David et al 



and the Centaurus cluster ( Tavlor et al.lll993 ). There is 
a correlation between the radio emission by the central 
galaxies and t he presence of a coo ling- fl ow that has been 
known since iBurns et al. (|l98lh and Ijones fc Forman 



(|l984D . The cooling gas is driven towards a massive ob- 
ject in the central galaxy that might produce the radio 
emission. However, what is quite surprising is that clusters 
exhibiting distinctive spiral-like structures are those with 
the greatest cooling-flows (see Table [2]) and containing a 
powerful radio ga laxy. The only exce ption is A2052 that 
contains an AGN ([Blanton et al.ll2003l . fjooih and does not 
exhibit a spiral-like substructure. 

The presence of a strong cooling-flow and a powerful 
central radio galaxy indicates that its radiation has not yet 
had time to heat the surrounding gas of the inner parts of 
the clusters. This can be understood if we consider a cycle 
in which the cool gas feeds the central AGN that begins 
to heat the surrounding gas. While the central gas is not 
sufficiently heated and the cooling flow is not suppressed, 
there will be a stage where both the AGN and the cool 
gas will coexist. 

The radio morphology of some clusters that exhibit 
spiral-like patterns appears quite different from most ra- 
dio galaxies be cause they fil l the X -ray cavities (regions 
of low density) . I Tavlor et al. |(|2Q02l) showed that the east- 
ern lobe of the central radio galaxy of Centaurus clus- 
ter appears to fill the region bounded by the spiral- 
like structure, while the western lobes turns around it. 
Similar confinement has been found for other radio galax- 
ies embended in dense ICM. For example, radio emission 
filling X-ray cavities was found at the core of Perseu s 
cluster teoehringer et al .111 9931 : iFabian et all 1 2 OOd [20021 ). 
A2029 (iTavlor et al.lll994h . and Hydra A cluster (A780) 



([McNamara et al.ll2000h . 



200 1[ ). Perseus ([Fabian et al.il2000l : IChurazov et al.ll2003h 



The ICM of merging clusters exhibit a variety of sig- 
natures of mergers such as cold- fronts and bow shocks. 
Although the most prominent cold fronts are produced in 
major mergers, cold fronts produced by minor mergers are 
expected to be far more frequent. Minor mergers generate 
subsonic turbulent gas motions within the cluster without 
destroying the cool-core. Different morphological features 
can be created by of minor me rgers, such as com et-like 
tails, bridges, plumes, and edges ([Poole et al.ll2006l ). Thus, 
one possibility for the origin of the spiral structure found 
in 7 clusters of our sample is that they are the signature 
of a small galaxy cluster or group minor merger with a 
more massive cluster. 

In a high-resolu tion numerical simulation of ideal- 
ized cluster mergers. lAscasibar fc MarkevitchI ( 2006[ ) con- 
cluded that if an infalling subcluster has a non-zero im- 
pact parameter, the cool gas of the main cluster acquires 
angular momentum, producing in a spiral-like feature of 
cold front due to the gas sloshing. In their simulation, af- 
ter 1.9 Gyr (see their Fig. 7) one can discern a spiral-like 
pattern in both the gas density and temperature distri- 
bution, that are morphologicaly similar to the structure 
found here. Since they used a mass ratio of 1:5, the cool- 
ing core is preserved. 
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The numerical simulations o fiG omez et al. also 
reproduce the results found here. In their merger 7, the 
subcluster has a small gas fraction and they assume a 
mass ratio of 1:16 so that the cooling- flow is not dis- 
rupted during the merger. An alternative explanation of 
the spiral-like structure is a contact discontinuity sepa- 
rating the main cluster gas from the gas of the infalling 
subcluster. However, it is also possible that the gas is 
stripped way from the infalling subcluster far from the 
core, never reaching t he inner region of the main cluster 
( Churazov et al.|[2QQ3l ). In this way, the observed strucutre 
consists only of the disturbed gas of the main cluster. A 
metallicity analysis would be able to distinguish between 
the above scenarios. If the cold gas in the spiral-like pat- 
tern is caused by the main cluster gas sloshing, the metal- 
licity in this region should be equal to the metallicity in 
the inner part of the main cluster. On the other hand, if 
it is a gas discontinuity, the metallicity in the spiral arm 
would be different from the inner part of the main clus- 
ter. The ICM enrichment analysis for this sample will be 
presented in a forthcoming work. 

In minor merger scenario, it is more likely that the 
cooling flow continues to feed the central galaxy. The ra- 
dio emission then, permeates the regions of low density, 
that is, fills the X-ray cavities because it is confined by 
the spiral-like feature. The ICM acts to distort the radio 
emission producing the radio morphologies, as can be seen 
in some clusters. 

We note that the detection of spiral-like features 
is not limited to the clusters in our sample but has 
also been achieved in other clusters, su ch a s A1975 

2004), A2029 
2009>). A1795 
1 l2009bh .This 



Fabian et al. l2QQlh 



Clarke et 



Jfl2004)) 



A262 (iBlanton et al 
A5098 (iRandall et Tl.! 



Liuzzo et a l. 2009), and A1201 (Q wers et a. 



kind of pattern has also been observed in g:a laxy groups 
(|Gastaldello et all l2007l : iRandall et al.ll2009h . indicating 
that they are not restricted to clusters. 

Since we find that a high number of clusters contain 
these structures, it would be interesting to bolster these 
findings by using numerical simulations. In particular, it 
would be instructive to study the frequency of off-axis 
merger halos, as a function of redshift, that produce a cool 
g:as spiral-like pattern when ad opting mass ratios of 1:5 
(lAscasibar fc MarkevitchI l2006l ) and 1:16 ( Gomez et al 



2002[ ). Unfortunately, these studies are beyond the scope 
of our present investigation. 

From the theoretical point of view, future numerical 
hydrodynamical simulations could detect the signatures 
of cluster mergers that could then be compared with ob- 
servations. In this sense, statistical studies of cluster mor- 
phology and substructures would provide important con- 
straints on models of structure formation and the evolu- 
tion of the observed cluster properties. Despite the large 
amount of available data, the dynamical evolution of ICM 
substructures remains poorly understood. Here we have 
studied in detail 15 nearby galaxy clusters highlighting 
the incidences of a particular ICM substructure morphol- 
ogy- 



6. Synthetic image 

We now present a synthetic image, created using the cfit- 
sio routine (a library of routines for reading and writing 
data files in the FITS data format). This synthetic im- 
age was constructed by adding multiple components to 
roughly mimic the Perseus cluster X-ray surface bright- 
ness. The synthetic image was the result of adding a 2D 
/3-model, an exponential spiral pattern, and two bubbles 
in the inner region of the cluster. The components and the 
result of their addition is seen in Fig. [3l 





Fig. 3: Synthetic image constructed to roughly reproduce 
the Perseus cluster X-ray surface brightness. Top panels 
from left to right: 2D /3-model, exponential spiral pat- 
tern, two bubbles, and the sum of these three compo- 
nents. Bottom panels from left to right: X-ray synthetic 
surface brightness, substructure map from this previous 
image, Perseus X-ray surface brightness, and its substruc- 
ture map. 



The goal of generating this synthetic image is to follow 
the same procedures applied to true clusters, on determin- 
ing substructure, and comparing with the substructure of 
the Perseus cluster. Although clearly present, one can see 
that the spiral pattern is not distinctive in the synthetic 
surface brightness image (see top panel of Fig. [3]), which 
shows the importance of this method in finding structures 
that are hardly discernible by eye. 

The synthetic-image substructure map shows the re- 
sult obtained in a spiral-like pattern structure very similar 
to the one found in the Perseus clusters, reinforcing the 
detectability of these spiral-like patterns. Furthermore, 
the result shows that the Perseus cluster X-ray surface 
brightness and substructures can be described well by 
the addition of three components, to produce in partic- 
ular an exponential spiral-like pattern, similar to that 
fo und in the results of numerical simulations performed 
bv lAscasibar MarkevitchI (|2006h . 



7. Conclusions 

In the most widely accepted scenario of structure forma- 
tion, clusters of galaxies form hierarchically by means of 
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mergers. Besides the modification of galaxy distribution 
and triggering of star formation, these merger events pro- 
duce variations in the physical properties, such as density 
and temperature variations that can be observed in X- 
rays. In this sense, the substructure maps have provided a 
powerful means of identifying hidden structures that are 
the last record of the merger history. 

We have presented 7 (A85, A426, A496, Hydra A clus- 
ter, Centaurus, Ophiuchus and A4059) out of 15 clusters 
that exhibit for spiral-like substructure correlated spa- 
tially with the boundaries of cold-fronts and cavities. The 
cavities inside these patterns are, in most of cases, associ- 
ated with radio emission. The aim of this study has been 
to determine the frequency of this morphological pattern 
at the centres of nearby clusters. Our analysis has shown 
how observational studies can improve the constraints on 
numerical simulations of structure formation and improve 
our understanding of ICM physical properties. 

Comparing our resu lts to numerical simulations 
( Ascasibar fc MarkevitchI l2QQ6l ). we have discussed the 
suitability of a model in which this spiral-like pattern is 
produced by an off-axis minor-merger that displaces the 
central cool gas. Since the minor merger does not disrupt 
the cool core, the central cool gas is still driven towards the 
central object that may produce the radio emission. This 
radio emission is bounded by the ICM inhomogeneities 
that will probably create an overall radio appearance that 
is quite different from most radio galaxies. Thus, we sug- 
gest that these distorted radio-emission morphologies may 
be caused by the confinement of dense X-ray gas in the 
spiral pattern. However, our interpretation of the spiral- 
like pattern and its relation to X-ray cavities and radio 
emission should be tested by additional more sophisticated 
simulations that ascertain possible merger signatures. 
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